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The extension of Majorana neutrino mass mechanism of the neutrinoless double-beta decay {Qufip) 
with the inclusion of right-handed leptonic and hadronic currents is revisited. While only the 
exchange of light neutrinos is assumed, the si /2 and pi/ 2 -states of emitted electrons as well as recoil 
corrections to the nucleon currents are taken into account. Within the standard approximations the 
decay rate is factorized into a sum of products of kinematical phase-space factors, nuclear matrix 
elements and the fundamental parameters that characterize the lepton number violation. Unlike in 
the previous treatments the induced pseudoscalar term of hadron current is included, resulting in 
additional nuclear matrix elements. An improved numerical computation of the phase-space factors 
is presented, based on the exact Dirac wave functions of the Si /2 and pi /2 electrons with finite nuclear 
size and electron screening taken into account. The dependence of values of these phase-space factors 
on the different approximation schemes used in evaluation of electron wave functions is discussed. 

The upper limits for effective neutrino mass and the parameters (A) and {rj) characterizing the 
right-handed current mechanism are deduced from data on the Oi^/3/3-decay of ^®Ge and ^®®Xe using 
nuclear matrix elements calculated within the nuclear shell model and quasiparticle random phase 
approximation. The differential decay rates, i.e. the angular correlations and the single electron 
energy distributions for various combinations of the total lepton number violating parameters that 
can help to disentangle the possible mechanism are described and discussed. 

PACS numbers: 23.40.-s, 21.60.Jz, 23.40.He 


I. INTRODUCTION 


The neutrinoless double-beta (0://3/3) decay is a process 
in which an atomic nucleus with Z protons decays to 
another one with two more protons and the same mass 
number A, by emitting two electrons and nothing else 

(A, Z) —i (A, -|- 2) -|- 2e . (1) 


Observing the Oiz/3/3-decay guaranties that neutrinos 
are massive Majorana particles - it means that neutrino is 
identical to its own antiparticle [T]. Oi/jS/J-decay violates 
total lepton number conservation and is forbidden in the 
Standard Model. 

When the light-neutrino exchange produced by left- 
handed currents is the driving mechanism for QvfJjS- 
decay, the relation between the effective Majorana neu¬ 
trino mass and the inverse half-life of the Ozz/3/3-decay can 
be written as [5] 


rr-iOlZ 
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( 2 ) 


where Z), and represent the phase-space 

factor, the axial-vector coupling constant and the nuclear 
matrix element of the process, respectively. In that case 
the ultimate goal of the search for Ojz/ 3/3-decay is the de¬ 
termination of the effective Majorana neutrino mass, 


TO/3/3 = C^e^TOi -b C/e2TO2 + C^esTOa. (3) 


Here, Uei and mi (i=l,2,3) are elements of Pontecorvo- 
Maki-Nakagawa-Sakata neutrino mixing matrix and 
masses of neutrinos, respectively. 

An improved calculation of G^^ in that case, taking 
into account the electron Dirac wave functions with finite 
nuclear size and electron screening, was performed in [3] . 
The main theoretical uncertainty is represented in the 
computed values of the nuclear matrix elements. There 
is a factor of 2-3 difference between the different methods 
of calculations of the . 

The left-right symmetric theories nig provide a nat¬ 
ural framework to understand the origin of neutrino Ma¬ 
jorana masses. In general one cannot predict the scale 
where the left-right symmetry is realized, but it is natu¬ 
ral to assume that it is as low as ^ a few TeV which can 
affect the Oi//3/3 decay rate significantly. 

In the left-right symmetric theories in addition to the 
left-handed V-A weak currents also leptonic and hadronic 
right-handed V-l-A weak currents are present. In that 
case new mechanism of the Ozz/3/3-decay need to be con¬ 
sidered. In the past the Oi//3/3-decay rate in the pres¬ 
ence of the right-handed leptonic and hadronic currents 
was discussed in pm. Recently, contributions to the 
Ojz/3/3-decay in a TeV-scale left-right symmetric models 
for type-I seesaw dominance were revisited [HHII]. By 
making a qualitative analysis without considering rel¬ 
evant phase-space factors and nuclear matrix elements 
it was found that Wl-Wr exchange (A mechanism) and 
Wr-Wr mixing (ry mechanism) could give dominant con- 
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tribution to the Oz^/3/3-decay amplitude by assuming a 
wide particle physics available parameter space includ¬ 
ing left-right neutrino mixing [iniiii]. We note that the 
discovery of a non-standard 0z//3/3-decay mechanism such 
as a right-handed current would rule out most models of 
baryogenesis at scales above 40 TeV [12]. 

The goal of this paper is to revisit the Oz^/3/3-decay 
mechanism due to the right-handed currents by consider¬ 
ing exact Dirac wave function of electrons and the higher 
order terms of nucleon current. We note that in |6] the ef¬ 
fect of induced pseudoscalar term of nucleon current was 
neglected and phase space factors were expressed using 
approximate electron wave functions for a uniform charge 
distribution in a nucleus by keeping only the lowest terms 
in the power expansion in r. In that context the subject 
of interest is the comparison of the power expansion ver¬ 
sus the exact treatment and the finite nuclear size effects. 
In this work the newly derived decay rate will be written 
in a compact form and the corresponding nuclear matrix 
elements will be presented by assuming the usual clo¬ 
sure approximation for intermediate nuclear states. The 
phase-space factors will be evaluated by using the exact 
Dirac wave functions with finite nuclear size and elec¬ 
tron screening. The differential characteristics, i.e. the 
angular correlations and the single electron energy dis¬ 
tributions will be described and discussed, and the decay 
rates and updated constraints on the lepton number vio¬ 
lating parameters for different combinations of the total 
lepton number violating parameters will be recalculated. 
This will make it possible to judge the importance of the 
A and t] mechanism. 


II. ELECTRON WAVE FUNCTIONS 


An important ingredient in the calculation of the elec¬ 
tron energy spectrum is the radial electron wave function 
distorted by the Coulomb field. We adopt the Dirac wave 
functions in a central field. 




9H.{e,r)xKt,i'r) \ 


(4) 


given, e.g. by Rose m- Here, e and r stand for energy 
and position vector of the electron, respectively, r = |r| 
and f = rjr. The index k takes positive and negative 
integer values (k = ±/c; k = 1,2,3-••). Total angular 
momentum is given as = |«;| —1/2 while orbital angular 
momentum takes values 


Ik = |k| - I 

= K 


for K < 0 

for At > 0 


(5) 


Radial wave functions gK{£,r) and /„(e,r) obey the ra¬ 
dial Dirac equations. 


dpK 

dr 


At -|- I 


Qk + {e-V{r) +me)fK 


XT = -(e - Lh) - H- 

dr 


where V (r) is the central Coulomb potential. The natural 
units h = c = 1 are used. 

The electron wave function expressed in terms of spher¬ 
ical waves is given by 

'k(e,r) = T’j-(-^ 


Here, superscript displays the orbital angular momen¬ 
tum (/„ = 0,1, 2, • • •) in a spectroscopic notation = 
s,p, d, •••) and the total angular momentum j„. The 
states of particular interest in our calculations are: 


V[/(si/2) 
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( 8 ) 


where p = p/p and p is the electron momentum. The 
asymptotic behavior of the radial wave functions for large 
values of pr is given by 


( gK{£,r) \ 

V fKi£,r) J 


(9) 


1 ( \l sin(pr - Z f -b 4 + aZf | log 2pr) 
\ \lcos{pr - -b 4 + log2pr) 


Zf is the charge of the final system which generates the 
potential V{r). 

In what follows, different approximation schemes for 
the calculation of radial wave functions g±i and f±i 
associated with emitted electron in the Si /2 and Pi /2 
wave states are briefly presented. 


The approximation scheme A): The relativistic elec¬ 
tron wave function in a uniform charge distribution in 
nucleus is considered. The lowest terms in the power 
expansion in r are taken into account. The radial wave 
functions take the form 


( g-ii.£,r) \ 

V J 


( g+i{£, 
V /-i(e> 




\jd^o{Zf,e) X 
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Here, Tfe_i (for fc = I, 2, • • •) is given by 


Fk-l = 


r(2fc-b I) 

r(fc)r(i + 27 fe) 
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where 


7/c = ^A;2 - {aZfY 

y = ( 12 ) 


This approximation scheme was commonly used in the 
past calculations of the phase-space integrals of double 
beta decay processes |S]. 

The approximation scheme B): The analytical solution 
of the Dirac equation for the point-like nucleus is consid¬ 
ered M- Radial wave functions take then the form 


5K(e,r-) = y — 


kI e + me\T{l + "/k+iy)\ 




= y — 


kpr\ 2e r(l-|- 27 fe) 

9 1^1 (7fe -iy,l + 27 fc,- 2 ipr)| , 

j e-m, |r(l + 7fc+»y)l (gp^p,.y /2 


kpr\ 2e r(l-|- 27 fe) 

sjj ^Fi{jk -iy,l + 27 fc,- 2 zpr)| , 

(13) 


with 


^ Ik-ty 


K — iyme/s 


(14) 


Here, iFi(a,h,z) is the confluent hypergeometric func¬ 
tion. 

The approximation scheme C): The exact Dirac wave 
functions with finite nuclear size, which is taken into ac¬ 
count by a uniform charge distribution in a sphere of 
nuclear radius i?, are assumed [3]. The numerical cal¬ 
culation can be accomplished by the subroutine package 
RADIAL [T^, where the input central potential is given 
by 


V{r) = 



for 

for 


r < R 
r > R 


(15) 


Here, R = with tq = 1.2 fm. In the code the ra¬ 

dial Dirac equations are solved by using piecewise exact 
power series expansion, which are summed up to a pre¬ 
scribed accuracy so that truncation errors can be avoided 
completely. 


The approximation scheme D): The exact Dirac wave 
functions with finite nuclear size and electron screening 
are used [3] . The effect of screening of atomic electrons is 
taken into account by the Thomas-Fermi approximation. 
It uses the solution of the Thomas-Fermi equation. 


(Pip y)3/2 

dx"^ yjx ’ 


(16) 



FIG. 1: The radial wave functions of an electron in the Si /2 
wave state, g-i{e) and /+i(£) (upper panels), and in the pi /2 
wave state, g+-i_{e), f-i{e) (lower panels), as function of the 
electron kinetic energy e — me. Results are presented for an 
electron emitted in the double /3-decay of ^®°Nd. In an evalu¬ 
ation of radial wave functions (w.f.) four different approxima¬ 
tion schemes are considered (see Sec.(|n| for details): A) The 
standard approximation of Doi et al! |6j; B) An analytical 
solution of Dirac equations for a pointlike nucleus is assumed; 
C) An exact solution of Dirac equations for a uniform charge 
distribution in nucleus is considered at nuclear surface; D) 
The same as the previous case but the electron screening is 
taken into account [3]. 


with X = r/h, where 


2/3 




( 17 ) 


The Thomas-Fermi function can be rewritten in terms 
of an effective charge p{x) = Zeff{x)jZf. Therefore, 
boundary conditions 


(^(0) = 1, 


Tioo) = — 


(18) 


of Eq. (161 take into account the fact that the final atom is 


a positive ion with electric charge -1-2. We adopt here the 
Majorana method described in Ref. m in order to solve 
Eq.(16). The input potential is then V{r) = p{r)Vo{r), 
where Vo(^) is defined in Eq. (15). 

In Fig. [2 radial wave functions of an electron in the 
Si /2 wave state {g-i{£) and /+i(e)) and in the pi /2 wave 
state ((/+i(e), /_i(e)) evaluated at r = ii are plotted 
as a function of the electron kinetic energy e — for 
the case of the double /3-decay of ^®°Nd. We see that 
wave functions A, which correspond to leading finite-size 
Coulomb, agree qualitatively well with wave functions of 
the other approaches in the case of (7-i(e) and /_i(e) 
but that there are significant differences for (/_|_i(e) and 
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f+i{e). These differences are apparent especially at lower 
energies. We notice also a rather good agreement be¬ 
tween results for wave functions of B and C in general. 
The screening of atomic electrons affects mostly the pi /2 
wave functions, but is essentially negligible for the Si /2 
states. 


III. DECAY RATE FOR THE NEUTRINOLESS 
DOUBLE-BETA DECAY 

One of the most prominent new physics model that 
incorporates the LNV and which leads to potentially ob¬ 
servable rates for the 0^/3/3-decay is the minimal left- 
right symmetric model (LRSM) [H which extends the 
standard model gauge symmetry to the group SU (2 )l ® 
SU{2)r 0 U{1)b-l. The right-handed neutrinos nec¬ 
essary appear here as a part of the SU{2)r doublets. 
The lepton multiplets Li = {vi,li) are characterized by 
the quantum numbers Qll = (1/2,0,—1) and Qlr = 
(0,l/2,-l) under SU{2)l ® SU{2)r ® U{1)b-l- The 
Higgs sector contains a bidoublet </> and two triplets 
and with vacuum expectation values (VEV) vr 
and vr, respectively. The VEVs fulhll the condition 
vlVr = . The VEV vr breaks SU(2)r 0 U{1)b-l 

to U{1)y and generates masses for the right-handed Wr 
and Zr gauge bosons, and the heavy neutrinos. The Wr 
and Wr are in general not mass eigenstates and are re¬ 
lated to the mass eigenstates Wi and W 2 with masses 
Ml and M 2 {Mi < M 2 ) as 

/ IV^ \ ^ / cosC sinC \ / Wf \ 

J [-sinC cosC J [W 2 - J- 

(19) 


Then, the effective current-current interaction which can 
trigger the Oiz/3/3-decay can be written as [6] 


ttP _ G/3 r . p , ■ p A 

^ ~ ^2 7 ^ TLp ^ a3l '^Rp 



+ VJrJrp + ^JrJrp + 


(20) 

Here, Gp = Gp cosOc, where Gp and 9c are Fermi 
stant and Cabbibo angle, respectively. The coupling 
stants A, 77 and x are chosen to be real. We have 

con- 

con- 

T] ~ -tanC, X = V, 

A ~ {Mwi/Mw2)'^- 


(21) 


]Vj, respectively. The electron neutrinos eigenstates can 
be expressed as 


3 

i'eL = ^ {Uejl'jR + Sej{NjR)^) , 

3=1 

3 

i^^R = ^ (t;. {,y,R)^ + v:^n,r) . (23) 

3=1 

The (3 -I- 3) scenario is assumed. The 3x3 block 
matrices in flavor space U, S, T, V, generalizations of 
the Pontecorvo-Maki-Nakagawa-Sakata matrix, consti¬ 
tute the 6x6 unitary neutrino mixing matrix m 

W = ( T y ) ^ (24) 


which diagonalizes the general 6x6 neutrino mass matrix 
in the basis (i/r, (Nr)^)"^: 


M 


Mr Mr, \ 
Ml Mr ) 


(25) 


with Majorana and Dirac mass terms, which are propor¬ 
tional to Mr » uhVl, Mr « UmVr and Mr, « yR,v, 
where yM and yo are the Yukawa couplings. The full 
parametrization of matrix U includes 15 rotational an¬ 
gles and 10 Dirac and 5 Majorana CP violating phases. 
It is possible to decompose U as follows m 


u 


{10 \ f A R \ f Vo 0\ 
{0 Uo J \ S B J {0 1 ) 


(26) 


where 0 and 1 are the 3x3 zero and identity matri¬ 
ces, respectively. The parametrization of matrices A, B, 
R and S and corresponding orthogonality relations are 
given in [T7]. In the limit case A = 1 , B = 1 , R = 0 
and S = 0 there would be a separate mixing of heavy 
and light neutrinos, which would participate only in left 
and right-handed currents, respectively. In that case only 
the neutrino mass mechanism of the Oj^/3/3-decay would 
be allowed and exchange neutrino momentum dependent 
mechanisms associated with the Wr-Wr exchange and 
Wr-Wr mixing would be forbidden. If masses of heavy 
neutrinos are above the TeV scale, the mixing angles 
responsible for mixing of light and heavy neutrinos are 
small. By neglecting the mixing between different genera¬ 
tions of light and heavy neutrinos A, B, R and S matrices 
can be approximated as follows: 


The left- and right-handed leptonic currents are given by 


Aw 1,R« 1,R« 

ruRRiv 


(27) 

turny 


Jl = e7p(l - lb)VeL, Jr = e7p(l -f l5)VeR. 

( 22 ) 

The VeL and UeR are the weak eigenstate electron neu¬ 
trinos, which are expressed as superpositions of the light 
and heavy mass eigenstate Majorana neutrinos Vj and 


Here, mR, represents energy scale of charge leptons and 
RT-LNY is file total lepton number violating scale, which 
corresponds to masses of heavy neutrinos. For sake of 
simplicity the same mixing angle is assumed for each 
generation of mixing of light and heavy neutrinos. We 
see that Uq can be identified to a good approximation 
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with the PMNS matrix and Vq is its analogue for heavy 
neutrino sector. Since Vq is unknown, it is common to 
assume that the structure of Vg is the same one as Uq- 
Assuming the non-relativistic impulse approximation 
the left and right hadronic currents and become 
0 


[{gv - gACn) 

n 

= [{gv + gACn) g^° 

71 

+ g’’^ i^g'A^n ~ g'v^n + g'p ^2rr^ ) 

(28) 


Here, q™ = Pn ~ Pn is the momentum transfer between 
the nucleons. The final proton (initial neutron) possesses 
energy E'^ (En) and momentum pj^ (p„). (T„, t+ and r„ 
are the Pauli matrix, the isospin raising operator and the 
position operator, respectively. These operators act on 
the n-th nucleon. 

The nucleon recoil operators Cn and D„ are given by 


a 


D, 


• (pn + p1) 

2m N 

(pn + P^) 

2m PI 



g ■ Qn 

2mjv ’ 


i (i + ^ X q» 

V gv ) 2TOAr 


(29) 


Here, mpp is the nucleon mass, qv = gv{q^), Qm = 
9m(9^), Qa = gA{q^) and qp = qp{q^) are, respec¬ 
tively, the vector, weak-magnetism, axial-vector and 
induced pseudoscalar form-factors in the case of left- 
handed hadronic currents. As the strong and electro¬ 
magnetic interactions conserves parity there are relations 
among form-factors entering the left-handed and right- 
handed hadronic currents [B]: 


gA ^g^ 
gv g'v’ 


5 m ^ 5 m g^ ^g^ 

gv g'v ’ gv g'v 


(30) 


We note that the induced pseudoscalar term of the space 
component of hadronic currents was not taken into ac¬ 
count in derivation of the Op/3/3-decay rate presented in 
[ 6 ]. This simplification is avoided here. 

Due to helicity matching of the propagating neutrino 
the decay amplitude can be divided into two parts: 

a) If both vertices are V — A ov V + A, the amplitude of 
the process is proportional to the neutrino mass mj. We 
shall denote the corresponding parts of the /3-decay 
amplitude L-L and R-R terms, respectively. 

b) If one vertex is H — A and the other is H -I- A, the 
four momentum of propagating neutrino q^ = (uj, q) con¬ 
tributes. The corresponding part of the amplitude, which 
is denoted as L-R, is further separated into two terms. 


the w-term and the q-term. 

In the case of L-L- and R-R-terms the dominant con¬ 
tribution is associated with the emission of electrons in 
the Si/ 2 -wave state [T5]. However, the q-term changes 
the parity and therefore it requires that one of the final 
electrons to be in si/ 2 -wave while the other to be in pi/ 2 - 
wave, or both electrons are in the Si/ 2 -wave and nucleon 
recoil operator is taken into account. Nevertheless, the 
q-term is not negligible since the w-term is suppressed by 
a factor £ 12/5 ~ 1/40 [ 6 ], that makes the q-term compa¬ 
rable or even larger in comparison with the w-term. 

The standard approximations of Ref. [B] are adopted: 

i) Only mechanisms with the exchange of light neutrinos 
are considered and contributions from heavier neutrinos 
are neglected. Recently, it was concluded in Ref. mm 
that mechanisms with the exchange of light neutrinos 
could give dominant contributions to the Qufip ampli¬ 
tude in a wide range of the LRSM parameter space. 

ii) Closure approximation. Within this approximation 
energies of intermediate nuclear states En — {Ei Ef) 12 
are replaced by an average of En — iEi + Ef)/2 ~ 10 MeV 
and the sum over intermediate states is taken by closure, 

iii) The R-R-pait of the amplitude, that is multiplied by 

factor |A^ becomes negligible in comparison 

with mjsp. Thus it is neglected. 

iv) The terms proportional to the square of the nucleon 
recoil operators are also neglected. 

v) For L-L-part of amplitude only electrons in the Si /2 
wave state are included. The inclusion of the pi /2 elec¬ 
trons leads only to negligible contribution to the Ov^f] 
standard decay rate [18]. 

vi) In the case of the L-R term, two-nucleon potentials 
associated with the spatial q and time ui components of 
neutrino exchange potentials are simplified as follows: 


c?q 


Ki^) = 




27t^ \q + A — S 12 + A + Ei 2 

»q.x 

q-{- A 

dq / I 


I 


£12 


27r2 \g-|-A-£i2 g-|-A-|-£i2 

dq 1 „i„.x 


{q + AY 




(31) 


where A = En — {Ei Ef)/2 and £12 = £1 — £ 2 - Here 
£1 and £2 represent the energies of the final electrons. 
Furthermore, contribution of the pi / 2 -wave electrons and 
term in which the nucleon recoil is multiplied by the small 
w-term are also neglected. 

vii) Since IxUeig'v/gv \ ^ Wejl, the coupling constant x 
in Hamiltonian (20 1 is neglected. 


viii) A factorization of phase-space factors and nuclear 
matrix elements is achieved by the approximation in 
which electron wave functions 5 ±i(£,r), f±i{£,r) are re¬ 
placed by their values at the nuclear radius R. The no- 
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tation 

5±i(e) = 6 '±i(e,^), /±i(e) =/±i(e, (32) 

is used. 


Within above approximations the Oz^/3/3-decay half-life 
can be written as 


rpQv 

^1/2 


-1 p 


Ou 


In 2 


= g\ \Mgt\^ I C„ 


\ me J 


+ Cm\ (A) cos V'l + Cmrj ( 77 ) COS V '2 

me me 

+ Caa (A)^ -h (r/)^ -h Cx^ (A) ( 77 ) COS {tpi - ^ 2 ) 


(33) 

The effective lepton number violating parameters associ¬ 
ated with right-handed currents and their relative phases 
are given by 


3 

(A) = A|^C/e,re*(gV/5v)|, 

3 

( 77 ) = g\Y,UejT:^l 

3 3 

7/>i = aig[(^mjU‘^j)C^UejT*j{g'y/gv)*], 
i=i i=i 

3 3 

7/>2 = arg[(^m,[/ 2 ,.)(^C/e,Te*)*]. 
i=i i=i 

(34) 


With help of (23) and by assuming (27), Uq ~ Vq and 
ig'vlgv) ^ 1 we get 


(A) « {MwjMw2 f ^ ICh 

mLNV 

( 77 ) « -tanC 1^1, (35) 

mLNV 


with 


ICI — IC23C12C13S13 C32C13 C13C23C32S13 

~Ci 2 Ci 3 {ci^Si 2 + S 33 ) I 

~ 0.82 (36) 

Here, Cy = cos(0y) and Sij = sin(0ij). ^ was evalu¬ 
ated by assuming the best fit values for mixing angles 
fi*i 2 , Oiz and 023 entering the PMNS matrix [19]. The ex¬ 
perimental upper bound on the mixing angle of left and 
right vector bosons is C < 0.013 and if the CP-violating 
phase in the mixing matrix for right-handed quarks are 
small one gets C, < 0.0025. The flavor and CP-violating 
processes of kaons and B-mesons make it possible to de¬ 
duce lower bound on the mass of the heavy vector boson 


Mw 2 > 2.9 TeV [TT|. In the LRSM there could be ad¬ 
ditional contributions to 07^/3/3-decay due to the double 
charged Higgs triplet. However, as pointed in Ref. m, 
in the considered case of type-I seesaw dominance, these 
contributions can be neglected. 

The coefficients Cj (I=mm, mA, mg, XX, 7777 and A 77 ) 
are expressed as combinations of nuclear matrix elements 
and phase-space factors: 

Cmm = (1 — XJ’+ Xt)^G'oi, 

C'mA = —(1 — XF + Xt) [X2-Go3 ~ Xl-I-G' 04 ] , 

CruT] = (1-XF+Xt) 

[X2+Go3 — X1-G'o4 — XpGoS + XF^oe] , 

C'aA = X 2-^02 + gXl + Goil — gXl + X 2 -GoiO; 

C'?7r7 = X 2 -|-G'o 2 + gXl-Goil — gXl-X2-|-GoiO + XpG'oS 
~ XpXrGqi + Xp^og, 

Cxr, = —2[x2-X2-|-Go2 — g (XI+X 2 + + X 2 -XI-) Goio 


+ gXl + Xl-Goil]. 


(37) 


The explicit form of nuclear matrix elements Mqt and 
their ratios xi is presented in section (HI B). The inte¬ 
grated kinematical factors are defined as 


Gok — 


Giml 


64775 In 2i?2 J Sie,+e 2 +Mf-M,) 

X {hok{£i,£ 2 ,R)cose + gok{ei,e 2 ,R)) 

X pip2£i£2d£ide2d{cos9) 

= rr^cos0+^' 


1-1 V in 2 


d{cos0). 


(38) 


where k = 1 , 2 , • • •, 11 . pi and P 2 are momenta of elec¬ 
trons and 6 is the angle between emitted electrons. The 
function s hokj ei, £ 2 , R) and goki£i,£ 2 , R) are defined in 
section ( HI Aj ). These definitions are independent of the 
weak axial-vector coupling constant g^- The quantities 
Gok are given in units of inverse years. We note that 
if the standard wave functions of electron (w.f. A) are 
assumed Ggio = G 03 and Gqii = G 04 . If in addition 
contributions from the induced pseudoscalar term of nu¬ 
cleon current are neglected the decay rate in Eq. ([3^ 


reduces to that given in [ 6 |. Quantity Gq^. is relevant for 
the angular correlation between the two electrons. We 


note that Gq 3 = Ggg = 0 . 


A. Components due to electron wave functions in 
the phase-space factors 

The Si /2 and pi /2 electron wave functions at nuclear 
surface associated with emission of both electrons enter 
into the phase-space factors through the functions pre¬ 
sented below. 
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For phase-space factors C?g^ related with the an¬ 
gular distribution of emitted electrons the quantities 
hQk{(^i,t 2 ,R) are: 

^01 = —4Csj;(ei)C'ss(e2), 

ho2 = -^C'ss(ei)C'ss(£2)) 

mi 

ho3 = 0, 

^04 = - 3 ^ ^ (C/p(^i)C'ss(e 2 ) + C'fp(£ 2 )C'ss(ei) 

+ Cf,(£2)aa(ei) + , 

^05 = (C'/p(^l)C'ss(e2) + C'ip(£2)C'ss(£l) 

+ Cfp(£2)as(£l) + Cfp(£l)C.«(e2)) , 

^06 = Oj 

“ 7-;;;2 + C'/p(£2)C'ss(£l) 

- Cfp(£2)a,(ei) - Cfp(£i)C.,(£2)) , 

/i08 = -^|^(C4(£i)Cfp(e2)+C4(£2)C'fp(£l) 

yTTl^TiJ 

+ C'ss(ei)C'pp(£2) + C'ss(£2)C'pp(£i)) , 

32 

^09 “ 7-;^C'ss(£i)C'ss(£2), 

(meR) 

how = --j^how — ho2, 

hoii = Qhoii + gho 2 + holo, (39) 

with 

how = 3 ^ 2 ^ (C/p(£l)Css(£ 2 ) - C'/p(£ 2 )C'ss(£i) 


“ 7-j;£ 2 (^i(£i)C'sp(£2) + C'+(£2)C' (£ 1 ) 

(Wei?) 

+ C'^s(£l)C'7'p(£2) + C'^s{£2)C'7'p(£i)) , 

^ 7- 77 ^ (“C'j)p(£l)C“ (£ 2 ) - C'pp(£ 2 )C'“(£i) 

(meR) 

+ C'j)p{£i)C'+(£2) -I-C'j)p(£2)C'+(£i) 

+ 2C-(£i)C-(£ 2 ) + 2C+(£i)C+p(£ 2 )) , 

509 = ^^(C+(£i)C'+(e2) + C-(£i)C-(£2)), 
(meR) 

9. 

5010 — — 2^010 — 502 , 

5011 = 95011 + g 502 + 5010 j (41) 

with 

5010 = (-C'+(£l)C'“ (£ 2 ) -I-C'+(£ 2 )C';),(£i) 

+ C'^s(£i)C' 7^(£2) - C'^s(£2)C'Jp(£i)) , 

5011 = 13^2^2 [C'pp(£l)C'ss(£2) + C'pp(£2)C'ss(£i) 

+ C'j)p(£i)C'+(£2) -I-C'j)p(£2)C'+(£i) 

- 2C-(£i)C'-(£ 2 ) + 2C+(£i)C+p(£ 2 )] . (42) 

Here, C,, are combinations of radial components of S 1/2 
and pif 2 wave functions, 

Css(e) = g-i(e)f+i(e), Cpp(si) = gi(e)f-i{e), 

Cip{e) = /_i(£)/+i(£), C'fp(£) = g_i(£)g+i(£), 

C±(£) = gi^(£)± fl^(e), C^p(£)=gl^(£)±fl^(e), 
Cfp{e) = g-i{£)f-i(£)±g+i{e)f+i{e)- (43) 


+ C%(e2)Ces(ei) - C'fp(£i)a«(£2)) , 

^011 = ( 3 ^ Jiy [*^4(^l)^/p(^2) + C'fp(£ 2 )C'fp(£i) 

+ C'ss(£i)C'pp(£ 2) -|- C'ss(£2)Cpp(£i)] . (40) 

In addition, the components goki^i, £ 2 , R) of the phase- 
space factors ([M|) are: 


B. Nuclear matrix elements entering the decay 
rate 


The expression for the 0z//3/3-decay half-life in (331 con¬ 
tains matrix element ratios xi £^od their linear combina¬ 
tions xi± s-od X 2 ±- The quantities xi £^re defined as 


501 — 5ii — C'i(£i)C'i(£ 2 ), 

502 = |^(C;(£l)Ci(£2)-C-(£i)C-(£2)), 

503 = — (C'i(£i)C-(£ 2 ) - Ci(£ 2 )C-(£ 1 )) , 

mg 

504 = 3^ Jl {~^ssi.£l)Csp{^‘2) - Cjs{^2)C~p(£i) 
+ C'+(£i)C' 7^(£2) -I-C'+(£2)C'7'p(£i)) - 503 / 9 , 

505 = (C'“ (£i)C';),(£2) -I-C” (£2)C^(£i) 

+ C'/7(£i)C' 7^(£2) -I-C'/7(£2)C'/p(£i)) , 

506 = (C'/;(ei)C'ss(£2) + C'+(£2)C',,g(£l)) , 


Xi = Mi/Mgt, (44) 

where I=F,T, wF, cuGT, wT, qF, qGT, qT, R and P 
and Mgt is the dominant Gamow-Teller matrix element 
associated with mechanism due to the left-handed 
currents. The combinations xi± and X 2 ± are given by 

Xl± = XqGT - QXqT ± ^XqF, 

X2± = XGTuj + XToj ± XFcj - gXlT- (45) 

The nuclear matrix elements Mj depend on the exchange 
potentials hi{r) through 

Mp^GT,T = ^ {Af \\hF^GT,T(r-)OF,GT,T\\ Ai) 




















fuF — fp, fuiGT — IgT 


fuiT — It- (49) 


Ma,F,uiGT,uiT = ^ {Af \\hi^F,uiGT,uiTir-)OF,GT,T\\ Ai) 


Mp = ^ilA 


hp{r_y+-+ 


_(r- X I-+) ^ 

^o ^r 


i?2 


AIqF,qGT,qT = ^ {Af |1/igF.gGT.qT (?'-)Of,GT,T || Ai) 

rs 

AIr = ^ {Af \\[hRG{r-)OGT + hRT{r-)OT]\\ Ai) , 


where Of^gt,t are the familiar operators l,i7i ■ (T 2 and 
3(cti • ri2){a2 ■ fi2). 

The two-nucleon exchange potentials hj{r) with = qF, 
qGT, qT, RG, RT and P can be written as 


where 


Here, r+ = (r^. -|- rs)/2, r_ = (r^ — rs). is the co¬ 
ordinate of decaying nucleon and ji{qr) (i=l,2,3) are the 
spherical Bessel functions. It is assumed that Pr+p'r — Oj 
Er — E) ~ Q and p,. — p). ~ q, where q is the momentum 
exchange. The form factors gv{q^)i 5a(9^), OAtiq^) and 
gv{q^) are defined in [20] and gA = 1.269. 

If right-handed currents are switched off, all terms in 
Eq.pSj) except that proportional to Ci vanish. The con¬ 
nection with the standard OF/3/3-decay formula (33) is 
then Goi = G^'' and Mgt(1 — Xf + Xt) = . 


hi{r) = ^ I //(g, 


q dq 

q + En — {Ei + Ef)j2 


(46) 


where 


Jgt = 


jo{g,r) 


9a 


fp = 


It = 


X (^5i(g^)- 

9W), 


9 A{q^)gp{q^) q^ , g|(g^) g" 


niN 


2 joiqr), 

9a 


4m^ 3 


j2{q,r) (gA{q^)gp{q^) q^ g%{q^)q‘^ 


9a 


K 


PIN 


4:mjq 3 


f 9vi9 ) . / ^ 

JqF = r - 2 -Jl(gr)g, 

9a 


, (3W)„ , pW) 9= 


-h 


9A{q^)gp{q^) q^ 


9\ 


tun 


rjiiqr), 


r r / /g\{q^) gp{q )gA{q ) 9 \ ■ f . 

fqT = - —q -- Ji{qr) 

3 VV 9a ^9a ruN 


9Uq^^ 


- 9 
fpG = 
fpT = 


2g\ 20to|, 


-R 
3m N 
-R 


6mN 


1-f 

1-f 


[2ji(gr)/3-j3(gr)] ) , 
9M{q'^)\ gA{q^)gv{q^) . 


9a 


9v{q^\ 

9M{q^ ) \ gA{q^)gv{q^ 
9v{q^ 


9a 


3o{qr)q^, 
j2{qr)q‘^, 


R^ gv{q^)gA{q^) . i ^ 
fp = -2-Ji(gr)g 




(47) 


And the two-nucleon exchange potentials hj{r) with 
1= wF, wGT and wT take the form 


= ~j 


g2 dq 


{q + E^-{E, + Ef)/2y 


(48) 


IV. PHASE-SPACE FACTORS WITH 
IMPROVED ACCURACY 

In Section |n| different ways of the treatment of radial 
wave functions g±i{e,r) and f±i{e,r) associated with 
emitted electrons in the Si /2 and pi /2 wave states were 
presented. The derivation of the Of (313-deca.y rate was ac¬ 
complished by considering electron wave functions for the 
point-like nucleus (wave function B) or an extended one 
(wave functions A, G and D); that allowed to separate 
phase-space factors and nuclear matrix elements. The 
accuracy of the calculation of phase space factors will be 
discussed next and the improved phase-space factors as¬ 
sociated with mechanisms due to right-handed currents 
obtained using screened exact finite-size Goulomb wave 
functions of emitted electrons (wave functions D) will be 
presented. 

A numerical computation of all II phase-space factors 
entering the OF/3/3-decay rate was performed by using 
previously described 4 types of wave functions (A, B, 
G and D) for a sample of 3 isotopes (^®Ge, ^^°Te and 
^®°Nd). Results are presented in Table We see that 
by using the standard treatment of electron wave func¬ 
tions, corresponding to leading finite-size Goulomb cor¬ 
rections (wave functions A), a significant difference with 
results of the other three approaches appear, especially 
for nuclei with large nuclear charge Z. Surprisingly, re¬ 
sults obtained with wave functions B, corresponding to 
an analytical solution of Dirac equations for a pointlike 
nucleus, better agree with results corresponding to wave 
functions G and D (exact solution of Dirac equations for 
a uniform nuclear charge distribution with the radius R) 
than those obtained by the standard treatment of wave 
functions (wave functions A). It indicates that the ex¬ 
act treatment of the Coulomb field plays more important 
role than the position of decaying nucleon in the nucleus. 
From the Table |l| it is apparent that effect of the screen¬ 
ing of atomic electrons on the wave functions of emitted 
electrons does not play an important role. 

The phase-factors differ among themselves significantly 
in magnitude. This fact is manifested in Fig. One can 
see that phase-factors obtained with standard wave func¬ 
tions (w.f. A) are always larger than those with phase 


































9 


TABLE I: Phase-space factors Gok (k=l, ■ • ■, 11) in units yr~^ for the 0z//?/3-decay of ^®Ge,^^°Te, and ^®®Nd. Calculation was 
performed by assuming different approximations for the radial wave functions g±i and f±i of an electron: A) The standard 
approximation of Doi et al. [^; B) An analytical solution of Dirac equations for a pointlike nucleus is assumed; C) An exact 
solution of Dirac equations for a uniform charge distribution in nucleus is considered; D) The same as the previous case but 
the electron screening is taken into account [^. 




"®Ge 




isoTe 




isoNd 



w.f. 

A 

B 

G 

D 

A 

B 

C 

D 

A 

B 

C 

D 

Goi.lQi^ 

0.261 

0.244 

0.240 

0.237 

1.807 

1.535 

1.453 

1.425 

8.827 

6.986 

6.432 

6.316 

Go2.101^ 

0.428 

0.404 

0.397 

0.391 

4.683 

4.064 

3.851 

3.761 

40.190 

32.401 

29.869 

29.187 

Goa.lO^® 

1.478 

1.340 

1.316 

1.305 

12.237 

9.566 

9.065 

8.967 

70.032 

49.465 

45.593 

45.130 

G 04 .IOI® 

0.501 

0.489 

0.477 

0.470 

3.625 

3.315 

3.086 

3.021 

18.343 

16.000 

14.348 

14.066 

Gos.lO^® 

0.791 

0.727 

0.572 

0.566 

6.390 

5.185 

3.842 

3.790 

28.537 

21.183 

15.061 

14.873 

Goe.lO^" 

0.605 

0.547 

0.536 

0.531 

3.091 

2.398 

2.258 

2.227 

11.922 

8.323 

7.591 

7.497 

Got.IO^® 

0.365 

0.345 

0.274 

0.270 

2.713 

2.383 

1.788 

1.755 

13.625 

11.362 

8.233 

8.085 

Go8.10“ 

0.245 

0.236 

0.151 

0.149 

2.877 

2.653 

1.579 

1.549 

16.833 

14.996 

8.564 

8.405 

Gog.lO^® 

1.360 

1.263 

1.238 

1.223 

6.398 

5.354 

5.063 

4.972 

27.582 

21.530 

19.799 

19.454 

Goio-lO^® 

1.478 

1.531 

1.423 

1.410 

12.237 

14.602 

11.616 

11.455 

70.032 

105.415 

72.249 

71.154 

Goii.lQi® 

0.501 

0.500 

0.484 

0.476 

3.625 

3.564 

3.220 

3.148 

18.343 

18.334 

15.376 

15.055 



FIG. 2: The phase factors Gok (k=l, 11) in units of 

yr~^ for the 0z^/?/3-decay of ^®®Nd. Results are presented 
for approximate electron wave functions (type A [ 6 ]) and ex¬ 
act Dirac wave functions with finite nuclear size and electron 
screening (type D [3]). The exponents p for k=l, ■ • 11 are 

14, 14, 15, 15, 13, 12, 10, 11, 10, 15, 15, respectively. 


factors calculated with the more advanced wave func¬ 
tions. Exact treatment reduces the value of all studied 
phase-space factors. 

The phase-space factors (see Eq. (38)) contain prod¬ 
ucts of g±i(e) and f±i{e), which were evaluated air = R 
from radial wave functions g±i{e,r) and f±i{s,r) (see 
Eq. (32)). Thus, it is assumed that the /3-decay of both 
nucleons happens at the nuclear surface. This procedure 
can be generalized with help of normalized to unity dis- 


TABLE II: The phase-space factor Gqi in units yr~^ for the 
0iz/3d-decay of ^®Ge, ^®®Te, and ^®®Xe. Results are presented 
for exact Dirac wave functions with finite nuclear size and 
electron screening at nuclear radius R (Exact, approximation 
scheme D) and those averaged over the distribution deduced 
from the analysis of the dominant nuclear matrix element 
(Exact and averaged - see Eqs. (@). 

Goi-10^^ 

76 Qe~ ^®°Te~~ ^®®Xe 
Exact 0.23681 1.42547 1.46187 

Exact and averaged 0.23987 1.47396 1.52851 


tribution function D{ri) as follows: 


roo 

g±i{e, R)= 9±i{£, ri)D{ri)dri, 

Jo 

pOO 

f±i{e,R) = / /±i(£,ri)H(ri)dri (50) 
Jo 


where in this particular case D{ri) = 6{ri — R). 

In Eig. the distribution function D{ri) (or equiv¬ 
alently D{r 2 )) is shown corresponding to the nuclear 
matrix element (associated with the mpp mecha¬ 

nism) of the 0i//3/3-decay of ^®Ge and ^^®Xe and calcu¬ 
lated within quasiparticle random phase approximation 
with restoration of isospin symmetry [2T]. We see that 
/3-decay of both nucleons happens mostly in the vicinity 
of nuclear surface. In Table [n| the phase space factor Gqi 
calculated with help of D{ri) = 5{rx — R) and D{ri) de¬ 
duced from calculated nuclear matrix elements of ^®Ge, 
^^°Te and ^^®Xe are compared. We see that the corre¬ 
sponding effect is very small for ^®Ge and is only about 
4-5 % in the case of ^^°Te and ^^®Xe. 
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FIG. 3: The normalized ri (r 2 ) dependence of for ^®Ge 
and ri and r 2 are absolute values of a position vector 

of /3-decaying nucleon in a nucleus. 


The improved phase-space factors Gqj (j=l, • • •, 11) in 
units yr~^ associated with left and right-handed mecha¬ 
nism of the 0i//3/3-decay for nuclei of experimental inter¬ 
est are presented in Table |III| They were obtained using 
screened exact finite-size Coulomb wave functions for Si /2 
and pi /2 states of electron (wave functions D). In Table 
the phase space factors Gqj associated with angular 


IV 


distribution of emitted electrons are presented. 


V. CONSTRAINTS ON THE EFFECTIVE 
TOTAL LEPTON NUMBER VIOLATING 
PARAMETERS 


Experimental Oiz/3/3-decay half-life limits may be used, 
in combination with the formula (33), to constrain the 
effective Majorana neutrino mass mpp and the effective 
coupling constants (A) and ( 77 ) of the right-handed cur¬ 
rents. This can be done provided the values of phase- 
space factors and nuclear matrix elements are available. 
We shall use the quasiparticle random phase approxima¬ 
tion (QRPA) [55] and interacting shell-model (ISM) [53] 
matrix elements for such analysis. Their values are pre¬ 
sented in Table |V] In the case of ISM the magnitude 
of matrix elements Mgt calculated in [53] are assumed. 
We note that these matrix elements were evaluated when 
the contribution from the induced pseudoscalar term of 
hadron current was not taken into account. In analy¬ 
sis below the case of CP conservation {ipi = '4’2 = 0) is 
assumed. 


Different contributions to the Oiz/3/3-decay rate (33) 
are associated with different products of effective lepton 


FIG. 4: The decomposition of coefficients Ci (I=mm, mA, 
mi], AA, rjr] and Xy, see Eqs. (331 and (371 ) on partial con¬ 
tributions associated with phase-spase factors Gok (k=l, 
• • -,11). The symbols standing for index I are shown on the 
x-axis. The partial contributions are identified by index k, 
which value is shown by the corresponding slope. Ratios 
Cj'^/Cl calculated with the ISM and QRPA matrix elements 
are displayed with left (in red) and right (in blue) slopes, re¬ 
spectively. Results for ^®Ge and ^®®Xe are presented in the 
lower b) and upper a) panels, respectively. 


number violating parameters TO/ 3 / 3 , (A) and ( 77 ), which 
values are unknown. The importance of these contribu¬ 
tions depends also on the value coefficients Gj (I=toto, 
toA, to? 7 , AA, 7777 and A 77 , which can be calculated. The 
quantity is a superposition of contributions Cj^ associ¬ 
ated with phase-spase factors Gok (k=l, • • -,11). In Fig. 
0we show ratios G'j^/Gj for the 0z//3/3-decay ^®Ge and 
^®Xe and both sets of nuclear matrix elements. We note 
that coefficients Cmm, G\\, Grjn, and GmTj are dominated 
by a single contribution associated with a different phase- 
factor. In the case of Gm\ and G\jf there is a competition 
of mostly two contributions. 


^and the 


III) we can 


Using these nuclear matrix elements (Table 
phase-space factors calculated here (see Table ^ 
deduce from the experimental data — ^-O x 10 ^® for 
’^®Ge decay [5S] and > 3.4 x 10^® for ^^®Xe decay 
[5i)] (we use here the combined limit from the EXO and 
KamLAND-Zen experiments) the constraints on the ef¬ 
fective right-handed current couplings (A), ( 77 ) and the 
effective Majorana neutrino mass mpp listed in Table 


VI The constraints in Table VI are of a similar mag¬ 
nitude as those in Table I of Ref. m- However, they 
are based now on the exact treatment of the phase-space 
factors as well as on the more complete account of nu¬ 
clear matrix elements. Fig. [^ shows the allowed regions 
for TO/ 3/3 and (A) (( 77 )) for ( 77 ) = 0 ((A) = 0). Re¬ 
sults are presented for the two sets of nuclear matrix 
elements (ISM [53] [53] and QRPA [55]) and the standard 
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TABLE III: Phase-space factors Goj (j = l, • • •, 11) in units yr~^ obtained using screened exact finite-size Coulomb wave 
functions for Si/2 and pi/2 states of electron (wave functions D). The Q values were taken from experiment when available, or 
from tables of recommended value [3]. Goi is associated with the mechanism generated just by mpis. In the case of dominance 
of the (A) ((??)) mechanism the decay rate includes phase factors G02, Goio and Goii (G02, Got, Gos, Goo, Goio and Goii). 
The remaining phase factors are due to interference of these mechanisms (see Eq. (|38[)). 



«Ga 

"6Ge 



““Mo 

llOpd 

““Cd 

“"Sn 

isoTe 

“®Xe 

“°Nd 

Qlpp [MeV] 4.27226 2.03904 2.99512 

3.35037 3.03440 

2.01785 

2.8135 

2.28697 2.52697 2.45783 3.37138 

Goi.lO^'* 

2.483 

0.237 

1.018 

2.062 

1.595 

0.483 

1.673 

0.906 

1.425 

1.462 

6.316 

Go2.10^'‘ 

16.229 

0.391 

3.529 

8.959 

5.787 

0.814 

5.349 

1.967 

3.761 

3.679 

29.187 

G03.IOI" 

18.907 

1.305 

6.913 

14.777 

10.974 

2.672 

11.128 

5.403 

8.967 

9.047 

45.130 

Go4.10^" 

5.327 

0.470 

2.141 

4.429 

3.400 

0.978 

3.569 

1.886 

3.021 

3.099 

14.066 

Go5.10^3 

3.007 

0.566 

2.004 

4.120 

3.484 

1.400 

4.060 

2.517 

3.790 

4.015 

14.873 

G06.IO12 

3.984 

0.531 

1.733 

3.043 

2.478 

0.934 

2.563 

1.543 

2.227 

2.275 

7.497 

Go7.101° 

2.682 

0.270 

1.163 

2.459 

1.927 

0.599 

2.062 

1.113 

1.755 

1.812 

8.085 

Gos.10^1 

1.109 

0.149 

0.708 

1.755 

1.420 

0.462 

1.703 

0.939 

1.549 

1.657 

8.405 

Go9.101° 

16.246 

1.223 

4.779 

8.619 

6.540 

1.939 

6.243 

3.301 

4.972 

4.956 

19.454 

Goio. 10^"^ 

2.116 

0.141 

0.801 

1.855 

1.359 

0.309 

1.418 

0.660 

1.146 

1.165 

7.115 

Goii.lO^® 

5.376 

0.476 

2.183 

4.557 

3.502 

1.010 

3.704 

1.955 

3.148 

3.238 

15.055 


TABLE IV: Phase space factors Gqj associated with angular distribution of emitted electrons (see Eq. |38l) in units yr~^ 
obtained using screened exact finite-size Coulomb wave functions for Si/2 and pi/2 states of electron (wave functions D). The 
Q-values of Tab. are assumed. 




«Ca 

"®Ge 

«“Se 

®®Zr 

“°Mo 

llOpd 

“®Cd 

““Sn 

laoTe 

“®Xe 

“°Nd 

■'^01 

.10“ 

8.010 

0.679 

3.141 

6.484 

4.951 

1.397 

5.153 

2.699 

4.328 

4.426 

20.101 

'-^02 

.10“ 

5.144 

0.113 

1.075 

2.769 

1.770 

0.239 

1.629 

0.587 

1.137 

1.111 

9.138 

-^04 

.10“ 

1.786 

0.152 

0.703 

1.456 

1.112 

0.314 

1.161 

0.608 

0.977 

1.000 

4.566 

'^05 

.10“ 

10.714 

0.910 

4.219 

8.734 

6.674 

1.884 

6.965 

3.650 

5.862 

6.002 

27.397 

p0 

.10“ 

8.458 

0.711 

3.422 

7.431 

5.706 

1.589 

6.026 

3.089 

5.016 

5.155 

24.824 

p0 

'-^08 

.10“ 

3.553 

0.402 

2.121 

5.383 

4.271 

1.251 

5.054 

2.651 

4.498 

4.787 

26.100 

p^ 

'^09 

.10“ 

5.024 

0.313 

1.379 

2.562 

1.904 

0.504 

1.795 

0.899 

1.397 

1.387 

5.899 

p0 

'-^010 

.10“ 

0.695 

0.028 

0.317 

1.118 

0.764 

0.114 

0.884 

0.334 

0.706 

0.741 

7.816 

p0 

■'^011 

.10“ 

1.790 

0.152 

0.707 

1.466 

1.120 

0.317 

1.172 

0.615 

0.988 

1.012 

4.594 


(w.f. A) and improved (w.f. D) description of electron 
wave functions. Note that limits on lepton number vi¬ 
olating parameters are softened a little when other lep¬ 
ton number violating parameters have non-vanishing val¬ 
ues at the same time in comparison with the case when 
only a single parameter is non-zero. By assuming C = 
0.013 and 0.0025 mentioned earlier and the current limit 
( 77 ) < 2.98 X 10“® (^^®Xe, ISM, w.f. D) we end up with 
iTiD/mLNv = 2.8 X 10“^ and 1.5 x 10“®, respectively. 
For Mw 2 = 2.9 TeV and (A) < 3.34 x 10"^ (^^SXe, ISM, 
w.f. D) we get mD/rriLMV = 5.0 x 10“®. Thus, from 
the more stringent limits on (77) we obtain 7?T,iArv/TeV 
= 0.3 - 2 m^i/MeV, in agreement with the assumption 
that the basic scale of LRSM is O(TeV). It is therefore 
obvious that already the present limits of Ozt/ 3/3-decay 


half-lives can be used to constrain meaningfully the al¬ 
lowed parameter space of LRSM, and that mechanism 
associated with right-handed currents can compete with 
the one based on m/ 3/3 that is so often used. 

VI. DIFFERENTIAL DECAY RATES FOR 
LIMITING CASES 

It is of interest to consider the angular correlations 
of the emitted electrons and the single electron energy 
spectrum for the three limiting cases of lepton number 
violating mechanism since with sufficient experimental 
accuracy one could distinguish between decays due to 
coupling to the left-handed and right-handed hadronic 
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TABLE V: Nuclear matrix elements and their ratios. The 
quasiparticle random phase approximation (QRPA) and in¬ 
teracting shell-model (ISM) matrix elements are from [22] and 
|23| . respectively. In the case of ISM matrix elements Mqt 
calculated in is used. 



76, 

Ge 

136 

Xe 

ISM 

QRPA 

ISM 

QRPA 

Mgt 

2.350 

3.014 

1.770 

1.120 

Xf 

-0.106 

-0.389 

-0.151 

-0.412 

Xi+ 

0.686 

0.811 

0.782 

1.969 

Xi- 

1.340 

2.917 

1.784 

4.052 

X2+ 

0.633 

0.302 

0.556 

0.229 

X2- 

0.912 

1.216 

0.965 

1.195 

XR 

0.684 

1.192 

0.955 

1.958 

Xp 

-0.544 

-0.176 

0.256 

-0.321 



FIG. 5: Limits on the effective neutrino mass mpis and right 
handed parameters rj (left panels, (A) = 0) and A (right pan¬ 
els, (rj) = 0) implied by the constraints on the Or//?/5-decay of 
^®Ge (lower panels, Ty 2 > 3.0 x 10^® [ 21 ]) and ^®®Xe (upper 
panels, Ty 2 — 3.4 x 10^® |26|L To derive the bounds, the val¬ 
ues of nuclear matrix elements calculated within the ISM |23| 
and the QRPA [22] are used. Results are presented for ap¬ 
proximate electron wave functions (type A) and exact Dirac 
wave functions with finite nuclear size and electron screen¬ 
ing (type D). Ellipsoids show the boundaries of the allowed 
domains. 

currents. It is assumed that some future Oiz/3/3-decay ex¬ 
periments, e.g. the SuperNEMO [27] or NEXT [28], will 
have a unique potential to measure the electron tracks 
and thus to observe the decay electron angular correla¬ 
tions and individual electron energy spectra. 


The differential rate for the 0+ —)■ O'*' 0z//3/3-decay with 
the energy of one of the emitted electrons £i ( £i is the 


TABLE VI: Upper bounds on the effective neutrino mass mpp 
and parameters {rj) and (A) associated with right-handed cur¬ 
rents mechanisms imposed by the constraints on the Ozz/5/5- 
decay of ^®Ge (upper panels, T^j 2 > 3.0 x 10^® [25]) and ^®®Xe 
{T 1 J 2 — 3.4 X 10^® [ 21 ]). Nuclear matrix elements of interact¬ 
ing shell-model (ISM) |23| {Mgt is from [ 21 ]) and quasipar¬ 
ticle random phase approximations (QRPA) [ 22 ] are used in 
analysis. CP conservation is assumed ('01 =02 = 0). The 
standard electron wave functions (w.f. A) [6] and screened 
exact finite-size Coulomb wave functions (w.f. D) are consid¬ 
ered. 








76, 

Ge 

136 

Xe 

w.f. 






A 

D 

A 

D 








QRPA 


\mpp\ 

[eV] 





0.321 

0.333 

0.285 

0.315 

\mpp\ 

[eV] 

(for ( 77 ) = 

= (A): 

= 0) 


0.271 

0.284 

0.251 

0.285 

iv) X 

10"® 





3.093 

3.239 

2.077 

2.337 

iv) X 

10"® 

[eV] (for 

iv) = 

:(A) = 

= 0) 

2.652 

2.807 

1.840 

2.118 

(A)x 

10-" 





4.943 

5.163 

3.822 

4.370 

(A)x 

10-" 

[eV] (for 

iv) = 

:(A) = 

= 0) 

4.841 

5.068 

3.792 

4.349 








ISM 


\mpp\ 

[eV] 





0.515 

0.535 

0.222 

0.245 

\mpp\ 

[eV] 

(for ( 77 ) = 

--{X): 

= 0) 


0.436 

0.458 

0.194 

0.220 

iv) X 

10"® 





6.370 

6.760 

2.975 

3.291 

iv) X 

10"® 

[eV] (for 

iv) = 

:(A) = 

= 0) 

5.464 

5.863 

2.628 

2.976 

(A)x 

10-" 





8.462 

8.841 

3.000 

3.378 

(A)x 

10-" 

[eV] (for 

iv) = 

:(A) = 

= 0) 

8.304 

8.694 

2.949 

3.336 


kinetic energy fraction with respect to Qpp of one elec¬ 
tron, i.e. £i = iiQpp -\- rrie and £2 = Q/3/3 + 2me — £ 1 ) 
and the angular distribution with the angle 6 between 
the two electrons for three limiting cases can be written 
as follows: 

i) Case mpp 0 ((A) = 0 and ( 77 ) = 0), 

dr = gWHcrf dCm-nr- (51) 

\ me J 

ii) Case (A) ^ 0 [mpp = 0 and ( 77 ) = 0), 

dV = g\\MGT? {X)"" dCxx- (52) 

hi) Case ( 77 ) 7 ^ 0 {mpp = 0 and (A) = 0), 

dV = g\\MGT? {v? dCr^r,, (53) 

where 

dCrnm ~ (1 Xf) ^^011 

1 2 

dC\x = X2-dGo2 + gXi+d^oii — gXi+X 2 -dGoio, 

1 2 

dCrjrj = X2+dGo2 + gXi-dGoii — gXi-X2+dGoio 






















13 




FIG. 6 : The single electron differential decay rate normalized 
to the total decay rate vs. the electron energy e (e = (e — 
'me)IQpp) for 0i//3/?-decay of ^®Ge (left panels a), c) and e)) 
and ^®®Xe (right panels b), d) and f)). Results are presented 
for \nipp\^ (panels a) and b)); (A) = ( 77 ) = 0), |(A)p (panels 
c) and d)); mpp = (rj) =0), 1 ( 77 )^ (panels e) and f)); mpp = 
(A) = 0) terms. An approximate electron wave functions (w.f. 
A) and exact Dirac wave functions with finite nuclear size and 
electron screening (w.f. D) are considered. Nuclear matrix 
elements calculated within the ISM |23| and the QRPA |22| 
are used in calculations. 


FIG. 7: The angular correlation factor (see Eq. (55l) vs. 
the electron energy e (e = (e — me)/Qpp) for Ot//?/ 3-decay 
of ^®Ge (left panels a), c) and e)) and ^^®Xe (right panels 
b), d) and f)). Results are presented for \mpp\'^ (panels a) 
and b)); (A) = ( 77 ) = 0), |(A)p (panels c) and d)); mpp = 
( 77 ) = 0), 1 ( 77 ) 1 ^ (panels e) and f)); mpp = (A) = 0) terms. An 
approximate electron wave functions (w.f. A) and exact Dirac 
wave functions with finite nuclear size and electron screening 
(w.f. D) are considered. Nuclear matrix elements calculated 
within the ISM [23] are used in calculations. 


+Xp'^^08 ~ XPXRdQoi + xjidGo9i 


Here oi/ao is the energy-dependent angular correlation 
coefficient, which depends also on the chosen limiting 
( 54 ) case for lepton number violating parameters. 


with 


dGok = 


X 


d cos ddii 


GjmlQpp 

647r^i?2 


{hok{£i,£ 2 ,R) cos 9 -I- gok{£i,S 2 , R))piP 2 £i £2 
aQ+aiCOs9 fc = 1, 2, • • •, II. 


(55) 


Here, Oq and are angular correlation coefficients. 
Gok = ln2 Gok- 

The differential decay rate can be written as 


dr 

dcosOdii 


In Fig. I^the single electron spectra normalized to the 
total decay rate are shown as function of the electron 
energy e for the 0p/?/3-decay of ^®Ge and ^^®Xe due to 
non-vanishing mpp, (A) and ( 77 ). This quantity, ideally 
accessible experimentally, depends only very weakly on 
the chosen isotope, set of calculated nuclear matrix ele¬ 
ments, and whether standard or improved description of 
electron wave functions is used. The different character¬ 
istics of these three limiting cases provide a possibility to 
identify which of the parameters is responsible for Op/ 3/3- 
decay. 


oo (1 -I- oi/oo cosd). 


(56) 
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FIG. 8 : The angular correlation factor (see Eq. ( |55| )) vs. the 
electron energy e (e = {e—m^jQfip) for Oi^/3/3-decay of 
Results are presented for |(A)|^ term [rn^^ = ( 77 ) = 0). An 
approximate electron wave functions (w.f. A) and exact Dirac 
wave functions with finite nuclear size and electron screening 
(w.f. D) are considered. Nuclear matrix elements calculated 
within the ISM [23] and the QRPA [22] are used in calcula¬ 
tions. 

In Fig. I^the angular correlation factors oi/oo are pre¬ 
sented as function of the electron energy e for the 0z//3/?- 
decay of ^®Ge and ^^®Xe due to non-vanishing TO/ 3 / 3 , (A) 
and ( 77 ). The ISM nuclear matrix elements are consid¬ 
ered. The results slightly depends on the type of elec¬ 
tron wave functions and manifest similar behavior for 
both isotopes. In Fig. the oi/ao behavior in detail 
for the (A) limiting case is shown. Note that results in 
this case are affected more significantly by the choice of 
nuclear matrix elements and that the choice of electron 
wave functions changes oi/oo only slightly for small and 
large values of electron energy. 

VII. SUMMARY AND CONCLUSIONS 

It is often assumed that if and when the neutrinoless 
double beta decay (Oj// 3/3) is observed, it will be caused by 
the exchange of a virtual light Majorana neutrino and its 
decay rate will be proportional to the square of the effec¬ 
tive Majorana mass It would be possible, therefore, 
to constrain or determine the magnitude of this funda¬ 
mental parameter based on experimental limits or values 
of the decay half-life. However, that is not the only possi¬ 
bility. Many other manifestations of the ‘Physics beyond 
the Standard Model’ that would cause the Oi//3/3 decay 
were considered in the past. Among them the possibility 
of the right-handed lepton and/or hadron currents that 
could perhaps compete with the mass mechanism was of¬ 
ten discussed, see e.g. mz]. Such possibility arises natu¬ 
rally in the left-right symmetric models. In that case the 


Ojz/ 3/3-decay half-life will depend not only on the mpp but 
also, perhaps dominantly, on the parameters that char¬ 
acterize the right-handed currents, denoted (A) and ( 77 ) 
here. 

When it is assumed that the right-handed currents ex¬ 
ist, the 0 i 7 / 3 ^-decay half-life can be expressed as a sum of 
products of the phase-space factors, nuclear matrix ele¬ 
ments, and the fundamental parameters that characterize 
the new physics. In this work the particular emphasis is 
on the reformulation of this relation, on a careful deriva¬ 
tion of all terms in that expression, and on the new and 
more general evaluation of the phase-space factors. The 
phase-space factors depend on the wave functions of the 
emitted electrons, and various approximations were used 
in the past in their calculation. We use here the ex¬ 
act solutions of the Dirac equation for the Si /2 and Pi /2 
electron states, solving the Dirac equation in the poten¬ 
tial that includes the nuclear finite size and the electron 
screening. The possible approximations to this problem 
are analyzed and discussed, and in particular it is shown 
that using just the first order expansion in r in order to 
include the nuclear finite size with the sufficient accuracy, 
is not enough. Complete table of accurate phase-space 
factors for nuclei of interest is given. Compared with 
the treatment that uses only the first terms in the r ex¬ 
pansion (denoted as approximation A in this work), the 
exact phase space factors (approximation D) are smaller, 
in particular in the heavier nuclei (^^°Te and ^®°Nd) the 
reduction is ^ 30 % or even more. 

It is also often assumed that the nucleons involved 
in the 07z/3/3-decay are essentially at the nuclear surface, 
hence the phase space factors are evaluated with the elec¬ 
trons placed at r = i?. The adequacy of that assumption 
was not tested until now. Here it is shown, see Fig. |^and 
Table [H] that it is a reasonable assumption, even though 
an increase in the phase-space factors by a few percent 
in the heavier candidate nuclei is expected. 

Having the full set of the phase-space factors, it is pos¬ 
sible by combining them with the full set of nuclear ma¬ 
trix elements evaluated elsewhere, to obtain simultaneous 
or separate limits for the fundamental parameters TO/ 3,3 
and those associated with the right-handed currents (A) 
and ( 77 ). It turns out that again the difference between 
the previously used approximation A (just the first term 
in the expansion in r) and the more exact treatment (ex¬ 
act Dirac electron wave functions with the nuclear radius 
R and electron screening) in the hnal limits is relatively 
benign in ^®Ce, enlarging the limits on the fundamental 
parameters only by about 5%. However, in the heavier 
nucleus ^^®Xe the effect is larger, 10-15%. 

It is well known that by convincingly determining the 
07Z/3/3 half-life one would obviously show that the total 
lepton number is not a conserved quantity. However, that 
determination by itself will be insufficient to decide which 
of the possible mechanism is responsible for the decay. If, 
in addition, the single electron energy spectra, and the 
angular distribution of the emitted electrons, could be 
detected, it will help substantially in that task. If one of 
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the possible parameters, (A), or ( 77 ) dominates, the 
single particle spectra and angular correlations will be 
a decisive tool to determine the mechanism. Formulas 
that determine these quantities and the corresponding 
phase space factors are shown here. In that case the exact 
treatment of nuclear size makes only little difference. 
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